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ABSTRACT 

v 

The  mechanical  Impedance  of  the  human  and  various  other  primate  species  heads 
was  determined  over  the  frequency  30-5000  hertz.  A  simple  model  was  developed 
that  closely  follows  the  observed  Impedance  characteristics.  Spring  and  damping 
constants  were  evaluated  and  comparisons  between  species  obtained.  An  Impact 
tolerance  curve  was  computed  based  on  the  model  predictions  with  a  maximum  strain 
criteria.  Various  Input  pulse  shapes  were  analyzed  and  the  effect  of  pulse  shape 
and  duration  studied  for  the  different  species.  Published  values  of  tolerable 
Impulses  were  examined  and  compared  with  the  model  predictions.  Some  discrepancies 
have  been  noted  and  analyzed.  Particular  emphasis  In  this  work  has  been  on  linear 
Impacts  to  the  side  of  the  head  where  angular  accelerations  do  not  predominate. 
Results  are  presented  in  the  form  of  tolerable  average  acceleration  versus  pulse 
duration  curves  fo*  various  pulse  shapes  derived  from  primate  head  Impacts.  A 
comparison  with  the  Wayne  State  tolerance  curve  and  Elband's  work  of  1959  has  also 
been  made. 
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INTRODUCTION 

A  critical  factor  In  the  design  of  a  motor  vehicle  Interior  Is  the  limiting 
deceleration-time  pulse  the  vehicle  can  undergo  without  occupant  Injury.  The 
type  of  tolerance  data  needed  for  this  Interior  design  Is  not  readily  available. 
Although  much  work  has  been  done  In  determining  tolerances  for  thoracic  and 
frontal  head  Impacts,  much  more  work  Is  needed  to  further  substantiate  existing 
data  and  to  generate  new  data  for  other  regions  of  the  body. 

The  acceleration-time  tolerance  curve  developed  at  Wayne  State  University 
over  the  past  several  decades  Is  the  only  empirically  based  tolerance  curve  avail 
able  for  frontal  head  Impacts.  Although  this  curve  Is  widely  used  as  a  standard 
for  head  Injury,  the  experimental  design  upon  which  It  Is  based  has  been  the  sub¬ 
ject  of  much  controversy.  The  original  Wayne  State  Tolerance  Curve  (W.S.T.) 
(Llssner  1960)  consists  of  six  data  points,  five  plotted  as  peak  acceleration  and 
one  arbitrarily  plotted  as  the  average  of  the  peak  acceleration  and  the  average 
acceleration. 

A  revised  W.S.T.  curve  was  Introduced  by  Patrick  (1963).  This  curve  Is  a 
composite  based  on  a  wide  variety  of  pulse  shpaes  and  striker  configurations. 

The  failure  criteria  used  was  generally  skull  fracture.  In  spite  of  the  many 
interpretive  difficulties  associated  with  this  curve.  It  has  been  In  the  past 
the  principal  source  of  hard  data  for  human  head  Injury  tolerance. 


The  work  upon  which  this  paper  Is  based  was  supported  In  part  by  a  contract  with 
the  Department  of  Transportation,  National  Highway  Safety  Bureau,  #FH-ll-7288  and 
a  contract  with  National  Institutes  of  Health,  National  Institutes  of  Neurological 
Diseases  and  Stroke,  #PH-43-67-1136. 
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Several  attempts  to  generalize  the  Wayne  State  Tolerance  Curve  have  been 
made  In  the  last  four  years.  One  such  characterization  of  this  curve  was  set 
forth  by  C.  W.  Gadd  (1966),  who  defines  the  severity  Index  as: 


S.I.  -  £T  an  dt  (1) 

where: 

a  ■  acceleration,  force,  or  pressure  of  response  function  producing  threshold 
or  Injury  of  given  degree, 
n  ■  weighing  factor  greater  than  1. 
t  ■  time,  seconds 

The  value  for  n  was  obtained  from  a  log- log  plot  of  the  Wayne  Curve  and 
Elband's  (1969)  work  for  splneward  accelerations  of  seated  humans.  The  limiting 
Index  for  survival  was  determined  from  points  on  these  curves.  This  severity 
Index  Is  tied  very  closely  to  pulse  shapes  In  setting  the  survival  limit. 

The  J  Tolerance  value  (Slattenschek  1968)  or  the  Vienna  Institute  Index  Is 
the  newest  generalization  of  the  Wayne  Curve.  In  this  effort,  a  s Ingle- degree- 
of- freedom  model  Is  assumed  and  the  model  constants  are  determined  by  fitting  a 
linear  differential  equation  to  the  Wayne  State  data.  A  J  Tolerance  value  Is  then 
defined  by: 


,  _  Xmax 
J  Xerta 


(2) 


where: 

Xerta  *  tolerable  amplitude  from  Wayne  State  Curve. 

Xmax  *  maximum  X  generated  by  putting  the  acceleration  pulse  In  question 
Into  the  model  differential  equation. 


thus: 


J  >  1  Not  survlvable 
J  <  1  Survlvable 
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Once  more  the  Index  Is  determined  from  the  Wayne  State  Tolerance  Curve. 

If  the  dynamic  response  of  the  head  Is  known  for  a  wide  range  of  frequencies 
and  If  Injury  can  be  related  to  this  response,  then  a  tolerance  curve  can  be  gen¬ 
erated  for  all  pulse  durations,  given  only  one  point  on  the  tolerance  curve, 
provided  the  Injury  mechanism  does  not  change  over  that  range.  This  paper  Is  an 
attempt  to  Investigate  this  hypothesis  from  a  model  based  on  driving  point  Imped¬ 
ance  data  from  human  and  other  primates. 

METHOD 

The  mechanical  driving  point  Impedance  of  human  and  other  primate  heads  was 
determined  over  a  frequency  range  of  30-5000  hertz  utilizing  the  following  exper¬ 
imental  design. 

The  primates  were  anesthetized  and  a  10  millimeter  (mm)  circular  hole  was 
cored  0.25  Inch  above  the  ear  canal  on  the  side  which  was  attached  to  an  electro¬ 
magnetic  shaker.  The  loading  fixture  was  then  fastened  to  the  skull  at  this  site. 
On  the  opposite  side  of  the  skull  a  similar  hole  was  made  and  a  miniature  accel¬ 
erometer  attached  (Figure  1). 

The  sub-human  primates  skull  via  the  load  cell  was  rigidly  attached  to  the 

platen  of  a  200-pound  electromagnetic  shaker  (Figure  2).  The  bodies  were  supported 

In  a  sling  hammock.  A  servo  controller  was  then  set  to  apply  a  sinusoidal  constant 

amplitude  acceleration  of  either  10  or  20  g's  to  the  head  (Figure  3).  In  addition, 

an  accelerometer  was  placed  on  the  free  side  of  the  head  and  the  transmitted  accel- 
(Stalnaker  and  McElhaney  1970) 

eratlon  recorded. +  A  sweep  oscillator  drove  the  shaker  system  over  a  30  to  5000  Hz 
cycle  range,  while  an  automatic  on-line  analogue  Impedance  computer  was  used  to 
convert  the  force-time  and  acceleration-time  Information  Into  a  phase  and  Impedance 
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Figure  2.  Mechanical  Impedance  Test  Set-Up 


versus  frequency  plot.  With  this  system,  the  test  could  be  performed  In  less 
than  one  minute  (depending  on  sweep  rate)  and  produces  a  continuous  plot  of 
Impedance.  Mechanical  Impedance  versus  frequency  on  an  x,  y1 ,  y2  recorder  was 
recorded  for  the  living  anesthetized  sub-human  primates  (Figure  4).  The  pri¬ 
mates  tested  were  then  sacrificed  and  the  test  procedure  repeated.  The  heads 
were  removed  from  the  bodies  and  the  same  parameters  recorded  again.  The  skin 
and  mandible  were  removed  and  the  test  repeated.  The  brains  were  then  removed 
through  the  foramen  magnum  and  the  mechanical  Impedance  and  acceleration  on  the 
free  side  of  the  skulls  were  recorded  (Stalnaker  and  McElhaney  1970). 

The  same  experiment  was  performed  on  an  unembalmed  71 -year  old  human  male 
cadaver,  who  had  been  dead  approximately  30  hours  prior  to  the  experiment.  Con¬ 
stant  accelerations  of  1,  5  and  10  g's  were  applied  over  the  frequency  range  30- 
5000  Hz  and  the  abov»  mentioned  mechanical  parameters  measured. 

The  two-degree-of- freedom  system  shown  In  Figure  5A  has  been  developed  to 
closely  approximate  the  Impedance  characteristics  of  the  heads  as  measured  In 
these  experiments. 

If  the  system  Is  represented  schematically  as  In  Figure  5B,  the  system 
elements  are  combined  In  parallel  and  series.  Using  the  rule  of  parallel  systems 
the  impedance  at  point  4  Is 

z4  -  z3  +  z2  (3) 

Using  the  rule  of  series  system  the  Impedance  z2  at  point  2  Is 


z2 


1 

Vzc 


(4) 
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Figure  4.  Mechanical  Impedance  of  Primate  Head 
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so  now 


(5) 


f + 


Substl tutlng 


or 


(6) 


(7) 


This  model  has  one  anti  resonance  and  one  resonance.  At  low  frequencies,  the  system 
Impedance  approximates  the  total  mass  of  the  system;  at  high  frequencies  It  approx¬ 
imates  the  Impedance  of  the  drive  mass  element  m^.  The  phase  angle  shifts  from  +90° 
through  0°  at  the  anti resonance  frequency  to  -90°,  and  from  -90°  through  0°  at 
resonant  frequency  back  to  +90°.  The  height  of  the  peak  and  the  depth  of  the  valley 
are  controlled  by  the  amount  of  damping.  The  spring  can  be  approximated  for  this 
model  by  stiffness  line  going  through  the  Inflexion  point  of  the  portion  of  the 
mechanical  Impedance  curve  between  the  antiresonance  and  resonance. 

This  theoretical  equation  was  programmed  for  the  1130  IBM  Digital  Computer. 

By  varying  the  magnitude  of  the  mass,  spring,  and  dash  pot  of  this  model  a  best 

Macaca  mulatta 

fit  of  the  mechanical  Impedance  between  the  model  and  the+test  data  can  be  approx¬ 
imated  (Figure  6).  (Stalnaker  and  McElhaney  1970) 

The  mechanical  Impedance  of  the  cadaver  head  was  similar  enough  to  that  of 
the  living  monkey  head  so  the  same  model  could  be  used.  The  model  constants  that 
provide  a  best  fit  with  the  cadaver  head  data  are  n^g  =  .4  lb,  m2g  =  9.0  lb,  c  * 

2.4  lb/1 n •  k  ■  2.6  x  104  lb-sec/1n.,  with  an  antiresonance  at  180  hertz  and  a 
resonance  at  820  hertz  (Stalnaker  and  McElhaney  1970). 
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In  Interpreting  the  model  constants,  the  following  considerations  apply. 

The  calvarium  Is  divided  Into  approximately  four  major  sections:  the  frontal 
bone,  left  and  right  parietal  bone,  and  the  occipital  bone,  which  In  the  case 
of  a  monkey  Is  almost  entirely  under  the  brain.  The  attachment  to  the  shaker 
was  made  through  one  of  the  parietal  sections.  These  sections  of  the  skull  are 
connected  by  sutures  which  provide  Isolation  from  one  section  to  another.  This 
Implies  that  m^  In  the  model  may  be  one  or  more  of  these  sections.  The  Macaca 
mulatta  parietal  sections  were  found  to  weigh  approximately  0.065  pounds,  thus 
It  Is  believed  that  m^  In  the  model  corresponds  approximately  to  the  parietal 
sections  of  the  skull.  Comparing  the  spring  constant  obtained  In  the  monkey 
series  Impedance  tests  (Stalnaker  and  McElhaney  1970)  (1.8  x  104  Ib/ln)  with  the 
spring  constant  obtained  through  the  Impedance  modeling  (104  Ib/ln)  leads  to  the 
conclusion  that  the  spring  element  In  the  model  corresponds  for  the  most  part  to 
the  skull  stiffness.  Comparing  the  damping  constant  for  the  whole  head  of  the 
Rhesus  monkey  (1.2  lb-sec/in)  with  that  of  the  skull  alone  (0.6  Ib-sec/ln)  Indi¬ 
cating  that  half  the  damping  Is  due  to  the  skin,  muscle  and  brain  (Figure  7). 

With  this  linear  two-degree-of- freedom  model  as  a  mathematical  analogy  of 
the  head,  many  dynamic  Inputs  to  the  head  can  be  studied.  The  model  response 
can  be  expressed  In  terms  of  the  following  linear  differential  equations  (Figure 
8): 


mjJtj  ■  c(*2  -  *j)  +  k(x2  -  Xj) 

(8) 

1112X2  *-c(*2  *  )  *  k(x2  *  X| ) 

(9) 
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Figure  7.  Model  of  the  Head 


where 


(10) 


X  ■  x2  -  Xj 


Afiu 


thus 


x2  ■  X  +  x1 


(10a) 


Letting 


mgX  +  cJ(  +  kX  ■  -mgXj 


(11) 


x1  -  a(t)  any  Input  acceleration 

then 

X  +  -^X  +  -^X»  -a(t)  (11a) 

m2  m2 

Substituting  equation  (10)  Into  equations  (8)  and  (9)  and  then  substituting 
equation  (9)  from  equation  (8)  we  finally  arrive  at 

X  +  (1  +  %  £  X  +  0  +  ^)  -  0  (12) 

m-j  m2  m^j  m2 

The  required  equations  of  motion  of  the  model  are  therefore  equation  (11a) 
for  a  forced  vibration  Input  and  equation  (12)  for  a  free  vibration.  With  these 
two  equations  and  the  model  constants  developed  above,  the  dynamic  response  of 
the  head  model  can  be  studied  for  a  variety  of  Input  Impulses. 

Human  tolerance  curves  have  been  developed  from  Impedance  data  for  lateral 
driving  point  Impedance  data  discussed  above  and  from  frontal  Impedance  data  com¬ 
piled  by  V.  R.  Hodgson  (1968).  In  addition,  the  primate  lateral  Impact  tolerance 
curves  were  obtained  from  the  lateral  Impedance  data.  In  this  way,  six  data  sets 
were  formed  (Table  I). 
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The  criterion  for  Injury  was  assumed  to  be  the  average  strain  In  the  brain. 
This  then  Is  the  X  In  the  model  normalized  by  the  linear  dimensions  of  the  brain 
at  the  point  of  impact  and  In  the  direction  of  the  Impact. 

The  maximum  deflection  X  for  human  heads  In  the  longitudinal  direction  was 
determined  from  A.  N.  Elband's  (1959)  work  on  abrupt  transverse  decelerations 
(Figure  9).  A  rectangular  pulse  of  50  G's  for  45  msec,  was  used  as  a  survival 
acceleration  pulse.  This  same  strain  used  for  the  longitudinal  direction  was 
assumed  as  the  limit  for  lateral  Impact,  the  underlying  assumption  here  Is  that 
the  brain  Is  equally  vulnerable  to  strain  In  all  directions. 

The  maximum  X  for  sub-human  primates  was  calculated  from  lateral  head  Impact 
data  (Figure  10).  The  sub-human  primates  were  Impacted  on  the  side  of  the  head 
with  a  rigid  constant  velocity  Impactor.  The  extent  of  the  Injuries  were  deter¬ 
mined  In  72-hour  post-impact  autopsy.  The  tolerable  Impact  level  for  each  primate 
was  determined.  The  acceleration- time  pulse  for  the  tolerable  Impact  was  then 
used  as  the  head  model  data  Input. 

RESULTS 

The  results  of  this  study  are  represented  In  the  form  of  average  acceleration 
versus  pulse  duration  curves. 

The  average  strain  generated  by  a  (50  G,  45  msec)  rectangular  pulse  In  the 
longitudinal  direction  In  the  human  head  was  found  to  be  0.00329  In/In.  For  the 
same  given  strain  level  a  15  g  spread  was  found  above  10  msec,  depending  on  the 
pulse  shape.  A  minimum  was  found  to  occur  with  a  tolerable  value  of  36  g's  at 
3.7  msec.  In  the  sine  pulse  Input,  and  a  value  of  30  g's  at  5.6  msec.  In  the  tri¬ 
angular  pulse  Input.  No  minimum  occurred  for  the  rectangular  pulse  curve  (Figure 
11). 
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Figure  10.  Squirrel  Monkey  Head  Impact  Set-Up 


Figure  11.  Maximum  Strain  Criterion  Curves  for  Human  Longitudinal  Head  Impact 


The  same  strain  level  was  used  to  set  the  tolerance  level  for  the  lateral 
Impacts.  The  shape  of  the  lateral  Impact  curves  were  found  to  be  similar  to 
the  longitudinal  Impact  curve  with  a  shift  of  the  time  and  acceleration  axis. 

The  spread  was  only  4  g's  above  15  msec.,  that  Is,  the  sine  flattens  off  at  26 
g's,  the  rectangular  at  25  g's  and  the  triangular  at  22  g's.  A  dip  of  17  g's 
at  5  msec,  occurred  for  the  sine  curve  and  a  dip  of  15  g's  at  7.1  msec,  for  the 
triangular,  with  no  dip  In  the  rectangular  curve  (Figure  12). 

The  sub-human  primate  tolerance  curves  were  developed  for  the  squirrel  mon¬ 
key,  and  the  Rhesus  monkey  from  Impact  tolerance  data.  Maximum  strain  levels  of 
0.089  In/In.  and  0.098  In/In.  respectively  were  obtained  (Figure  13). 

CONCLUSIONS 

The  conclusions  of  this  study  Indicates  that  there  ‘s  a  significant  differ¬ 
ence  between  the  Wayne  State  Tolerance  Curve  and  the  Maxlmjm  Strain  Criterion 
(MSC)  presented  here. 

The  shape  of  the  input  pulse  significantly  effects  the  tolerable  accelera¬ 
tion  for  durations  greater  than  1  msec. 

Since  av  Is  dependent  on  the  pulse  duration  and  shape,  a  sine  pulse  Is  more 
"tolerable",  for  large  pulse  duration,  while  a  rectangular  pulse  Is  "tolerable" 
for  short  pulse  duration. 

For  pulse  durations  less  than  1  msec  the  shape  and  location  of  the  MSC  curves 
are  Independent  only  on  the  model  constants. 

The  human  lateral  MSC  curves  were  found  to  be  50%  lower  than  the  human  longi¬ 
tudinal  MSC  curve. 

The  squirrel  monkey  lateral  MSC  curve  has  the  highest  tolerance  level  with 
the  Rhesus  monkey  next.  The  tolerable  strain  level  for  both  animals  were  found 
to  be  quite  similar.  Thus  atleast  for  these  two  sub-human  primates  the  MSC 
tolerance  concept  holds. 
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Figure  12.  Maximum  Strain  Criterion  Curves  for  Human  Lateral  Head  Impact 
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NOMENCLATURE 


a  -  acceleration,  ft/sec2 
c  *  damping  constant,  lb-sec/ln. 

F  *  force,  lb 

i-fT 

k  *  spring  constnat,  lb/in. 
m  -  mass,  slugs 
t  *  time,  sec. 
v  *  velocity,  Ips 
x  ■  displacement.  In. 
z  »  mechanical  Impedance,  lb-sec/ln. 
e  -  strain  rate,  1/sec 
a  *  stress,  Ib/sq  In. 

<t>  *  phase  angle,  deg. 
u  *  frequency,  rad/sec 
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